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SYSTEM AND METHOD FOR DETERMINING A FLOW PROFILE IN A 
DEVIATED INJECTION WELL 

CROSS-REFERENCE TO RELATED APPLICATIONS 

This application claims priority from US Provisional Application 60/510,596, 
filed October 10, 2003, which is incorporated herein by reference. 

BACKGROUND OF THE INVENTION 

Field of the Invention 

[0001] The present invention relates to a system and methodology for 

determining a flow profile in a well, and particularly to determining a flow profile in a 
deviated injection well. 

Description of Related Art 

[0002] In a variety of wells, various parameters are measured to determine 

specific well characteristics. For example, temperature logging has been used for 
profiling the injection rate in vertical wells. Existing methods of analyzing injection 
profiles are designed for vertical wells where the injection interval is usually small and 
the time to flush the wellbore volume is negligible. Also, the displacement process of the 
reservoir fluid can be represented by a radial flow model. 

[0003] However, if the wellbore is deviated, such methods do not enable profiling 

of the injection rate. Thus, deviated wellbores, such as horizontal wellbores, present 
greater problems in evaluating and predicting flow profiles for injection wells. 
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BRIEF SUMMARY OF THE INVENTION 

[0004] In general, the present invention provides a system and methodology for 

using a well model in determining characteristics of an injection well. For example, the 
system and methodology enables the use of temperature profiles in a deviated injection 
well to determine a flow profile in such well. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0005] 1 Certain embodiments of the invention will hereafter be described with 

reference to the accompanying drawings, wherein like reference numerals denote like 
elements, and: 

[0006] Figure 1 is an elevation view of a completion and sensing system deployed 

in a deviated wellbore, according to an embodiment of the present invention; 

[0007] Figure 2 is an elevation view of the system illustrated in Figure 1 

illustrating a fluid being injected into the well, according to an embodiment of the present 
invention; 

[0008] Figure 3 is a flowchart generally representing an embodiment of the 

methodology used in determining a flow profile in a well, according to an embodiment of 
the present invention; 

[0009] Figure 4 is a diagrammatic representation of a processor-based control 

system that can be used to carry out all or part of the methodology for determining flow 
profile in a given well, according to an embodiment of the present invention; 
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[0010] Figure 5 is a flowchart generally representing a methodology for 

detennining flow profiles based on temperature profiles during fluid injection into a 
deviated well, according to an embodiment of the present invention; 

[0011] Figure 6 is a graphical representation plotting temperature against distance 

along a wellbore during an early period of injection; 

[0012] Figure 7 is a graphical representation similar to that of Figure 6 but with a 

different injection geometry; 

[0013] Figure 8 is a graphical representation similar to that of Figure 6 but at a 

later injection time; 

[0014] Figure 9 is a graphical representation similar to that of Figure 7 but at a 

later injection time; 

[0015] Figure 10 is a flowchart generally representing a methodology for 

determining flow profiles based on temperature profiles in a deviated well during a shut- 
in period, according to an embodiment of the present invention; 

[0016] Figure 1 1 is a graphical representation plotting temperature against 

distance along a wellbore during a shut-in period; 

[0017] Figure 12 is a graphical representation similar to that of Figure 1 1 but with 

a different injection geometry; 

[0018] Figure 13 is a schematic representation of a processor system that receives 

data related to temperature profiles and other well parameters to derive flow profiles; 

[0019] Figure 14 is a schematic representation of a deviated well divided into a 

multi-segment grid system for modeling; 
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[0020] Figure 15 is a graphical representation of dimensionless temperature 

plotted against dimensionless time; and 

[0021] Figure 16 is a flowchart generally representing a methodology for 

determining flow rates for a plurality of intervals along a deviated well. 

DETAILED DESCRIPTION OF THE INVENTION 

[0022] In the following description, numerous details are set forth to provide an 

understanding of the present invention. However, it will be understood by those of 
ordinary skill in the art that the present invention may be practiced without these details 
and that numerous variations or modifications from the described embodiments may be 
possible. 

[0023] The present invention generally relates to a system and method for 

determining flow profiles in a deviated well. Temperature measurements are taken along 
a wellbore, and those measurements are used in determining flow profiles along a 
deviated injection well, such as a generally horizontal injection well. In some 
applications, a flow profile is derived based on data obtained during injection of a fluid 
into the deviated well. In other applications, a flow profile is derived based on data 
obtained during a shut-in period following injection or during periods of resumed 
injection. 

[0024] A temperature sensing system, such as a distributed temperature sensor, is 

deployed with an operational completion and enables temperature measurements to be 
taken during fluid injection periods or during shut-in periods. Based on the collected 
temperature data, flow profiles of the injected fluid along the deviated well can be 
derived. 
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[0025] In general, when a cool fluid, such as a liquid, e.g. water or oil, or a gas, is 

injected into a hot reservoir, a variety of thermal changes occur. For example, during 
injection, cool fluid moves through the wellbore and into the reservoir while heat flows 
from the reservoir toward the wellbore. A similar effect occurs along the axis of the 
wellbore as fluid flows from the heel of the wellbore toward the toe, and heat flows from 
the toe of the wellbore toward the heel. The thermal characteristics of the heat flow can 
be modeled in a maimer that enables determination of the flow profile of the fluid flow 
into the reservoir. Other factors, such as thermal conductivity of the surrounding 
formation, may also be utilized in modeling the flow profile, as discussed below. 

[0026] Furthermore, if the injection of cool fluid is stopped, thereby creating a 

shut-in period, other unique thermal changes occur that enable the determination of 
injection profiles. For example, when the injection stops, the wellbore begins heating up, 
but not necessarily uniformly. The temperature recovery provides an indication of where 
the cool fluid was moving into the reservoir during injection. Reservoir intervals that 
were receiving substantial amounts of injection fluid, for example, are slower to rise in 
temperature during the shut-in period. These thermal changes are applied to models that 
enable the derivation of flow profiles, as discussed below. 

[0027] Referring generally to Figure 1, a system 20 is illustrated in accordance 

with an embodiment of the present invention. System 20 comprises a completion 22 
deployed in a well 24. In this example, well 24 is a deviated well having a generally 
vertical section 26 and a deviated section 28, such as the generally horizontal section 
illustrated in Figure 1 . Well 24 is defined by a wellbore 30 drilled in a formation 32 
having, for example, one or more fluids, such as oil and water. A tubing 33 extends 
downwardly into wellbore 30 from a wellhead 34 disposed, for example, along a seabed 
floor or a surface of the earth 36. In the illustrated example, tubing 33 extends to a casing 
shoe 38 that may be located at the lower end of vertical section 26 above a heel 40 of 
deviated section 28. Completion 22 is disposed in deviated section 28 and may extend 
from casing shoe 38 through heel 40 toward a toe 42 of well 24. In many applications, 



5 



WO 2005/035944 



PCT/IB2004/002877 



wellbore 30 is lined with a casing 44 that may be perforated to enable fluid flow 
therethrough. 

[0028] As further illustrated, system 20 comprises a temperature sensing system 

46. For example, temperature sensing system 46 may comprise a distributed temperature 
sensor (DTS) 48 that is able to continually sense temperature along deviated section 28 of 
wellbore 30 at multiple locations. Distributed temperature sensor 48 may be coupled to a 
controller 50 able to receive and process the temperature data obtained along wellbore 30. 
As discussed in greater detail below, controller 50 also enables use of the temperature 
data in conjunction with a model of the well to derive injection flow profiles of fluid 
flowing from completion 22 into formation 32 along deviated section 28 of the well 24. 

[0029] During injection, a fluid, such as water, is pumped down through tubing 

33 and into completion 22 along the illustrated horizontal section of the well. The fluid is 
forced outwardly along deviated section 28 such that fluid flows from completion 22 into 
formation 32, as indicated by arrows 52 of Figure 2. The injectivity profile or flow 
profile of fluid moving from the wellbore 30 into formation 32 often is not uniform. For 
example, the flow may be substantially greater in proximity to heel 40 as opposed to toe 
42. Also, the formation material may vary along the length of deviated section 28, which 
can extend a substantial distance, e.g. up to several kilometers. 

[0030] Referring generally to Figure 3, an example of the methodology of the 

present invention is illustrated in flow chart form. Determining flow profiles along a 
given injection well comprises deploying a sensor system in the well with an operable 
injection completion, as illustrated by block 54. Then, an injection fluid is injected into 
formation 32 via completion 22, as illustrated by block 56. The sensor system may 
comprise a distributed temperature sensor designed to sense temperature along deviated 
wellbore section 28, as illustrated by block 58. As discussed more fully below, the 
sensing of well parameters can be done during injection and/or subsequent to injection 
during a shut-in period. A well model may then be applied to determine the flow profile 
of injected fluid along deviated section 28 of well 24, as illustrated by block 60. 
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[0031] Some or all of the methodology outlined with reference to Figures 1-3 

may be carried out by controller 50 which comprises an automated system 62, such as the 
processing system diagramatically illustrated in Figure 4. Automated system 62 may be a 
computer-based system having a central processing unit (CPU) 64. CPU 64 may be 
operatively coupled to temperature sensing system 46, a memory 66, an input device 68, 
and an output device 70. Input device 68 may comprise a variety of devices, such as a 
keyboard, mouse, voice-recognition unit, touchscreen, other input devices, or 
combinations of such devices. Output device 70 may comprise a visual and/or audio 
output device, such as a monitor having a graphical user interface. Additionally, the 
processing may be done on a single device or multiple devices at the well location, away 
from the well location, or with some devices located at the well and other devices located 
remotely. 

[0032] In automatically determining a flow profile or profiles of fluid injected 

into formation 32, a model utilizing temperature changes along deviated section 28 as an 
indicator of flow profiles may be stored by automated system 62 in, for example, memory 
66. As illustrated best in Figure 5, the general approach involves obtaining an initial 
temperature profile along at least deviated well section 28, as indicated by block 72. 
Additionally, the temperature of the injected fluid, e.g. water, is measured, as indicated 
by block 74. This value also may be stored on automated system 62 for use in modeling 
the injectivity profile. In many applications of the modeling technique, a greater contrast 
between the temperature of the injected fluid and the temperature of the reservoir can 
improve the usefulness of the model. By way of example, the injected fluid may be at a 
temperature of 60-70 degrees Fahrenheit, and a reservoir may be at a temperature of 200- 
240 degrees Fahrenheit. However, the temperature of the injected fluid and the reservoir 
can vary substantially from one application and environment to another. The fluid is then 
injected, as illustrated by block 76. Subsequently, temperature profiles are obtained 
along the deviated section 28, and this data is provided to automated system 62 via the 
temperature sensing system 46, as illustrated by block 78. The temperature profiles can 
be taken during injection or during a shut-in period subsequent to injection depending on 
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the particular model applied. In either event, the model is used to derive a flow profile of 
the injected fluid, as illustrated by block 80. For example, the data collected may be 
processed according to the model/algorithm stored on automated system 62 to 
automatically present a well operator with detailed information on the injection flow 
profile via, for example, output device 70. 

[0033J In determining the flow profiles, a given a well model may utilize thermal 

behavior characteristics that occur during injection. An example of a well in which 
injectivity decays along a horizontal well axis is illustrated in Figure 6 which provides a 
graph of temperature plotted against distance along the horizontal section of the wellbore. 
The graph illustrates temperature changes along the horizontal wellbore section during 
the first few hours of injection. As illustrated by a graph line 82 at 0.002 day, graph line 
84 at 0.01 day, graph line 86 at 0.03 day, graph line 88 at 0.05 day, graph line 90 at 0.07 
day, graph line 92 at 0.1 day, and graph line 94 at 0.2 day, a gradual propagation of the 
temperature front in the wellbore is observed due to a much greater injectivity towards 
heel 40 than toe 42* A cooling or reduction of the temperature at heel 40 is also 
graphically illustrated. 

[0034] The thermal characteristics of a well in which injectivity is skewed 

towards toe 42 is illustrated in Figure 7. In this example a much faster propagation of the 
wellbore temperature front is observed, as indicated by the location of graph lines 82, 84, 
86, 88, 90, 92, and 94. 

[0035] Referring generally to Figures 8 and 9, temperature profiles for the 

examples provided in Figures 6 and 7, respectively, are again graphically illustrated, but 
at a substantially later period of the injection. In each of these figures, a graph line 96 
represents a temperature profile after one day of injection, and a graph line 98 represents 
a temperature profile after two days of injection. The thermal characteristics, such as 
those illustrated in Figures 6, 7, 8 and 9, demonstrate thermal changes that occur during 
injection of fluid into formation 32. Those thermal changes can be used by an 
appropriate model to determine flow profiles. In other words, the injection of a cooler 
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fluid into the reservoir at different rates along deviated section 28 creates thermal 
changes over the period of injection. By accurately measuring the thermal changes via, 
for example, distributed temperature sensor 48, the actual injection flow profiles can be 
derived by the appropriate model. 

[0036] At other phases of the process, useful thermal data also may be obtained. 

For example, following an injection phase, a shut-in phase can lead to interesting thermal 
events which can be modeled to provide an injection flow profile. When the injection of 
fluid starts, the wellbore begins heating, but not necessarily uniformly. The taking of 
temperature profiles during this temperature recovery period provides an indication of 
where the cooler injection fluid is moving into the reservoir during injection. For 
example, reservoir intervals receiving a greater flow of the cooler fluid are slower to 
regain heat during the shut-in period. Thus, the temperature profiles taken during a shut- 
in period can be used to determine injection flow profiles. 

[0037] As illustrated best in Figure 10, the general methodology for utilizing 

shut-in data involves obtaining an initial temperature profile along at least deviated well 
section 28, as indicated by block 100. The temperature of the injected fluid, e.g. water, 
also is measured, as indicated by block 102. The injection fluid is then injected, as 
illustrated by block 104. Following an injection period, e.g. two days of injection, the 
injection is shut down for a shut-in period, as illustrated by block 106. Subsequently, 
temperature profiles are obtained along the deviated section 28, as illustrated by block 
108. This data, along with other collected data, may be provided to automated system 62 
via the temperature sensing system 46. By applying the appropriate model, the shut-in 
temperature data can be utilized in deriving an injection flow profile along the deviated 
section 28 of well 24. In some applications, the injection is resumed, and that resumed 
injection may be followed by a subsequent shut-in period, as indicated by block 110. The 
repeating of injection and shut-in periods can be used to obtain additional data, to verify 
results, and/or to continually monitor the injection flow profile. It should also be noted 
that appropriate models can also be designed to utilize the thermal characteristics of a 
well when injection is resumed after a shut-in period. As the fluid is reinjected, the well 
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warms up and a step rise in temperature is indicated. This step or slug moves as a front 
along the deviated section of the wellbore and provides an indication of the flow profile. 
If; for example, the front moves slowly, this generally indicates greater flow towards the 
heel of the deviated section. If, on the other hand, the front moves more rapidly, this can 
indicate greater flow toward the toe of the wellbore. 

[0038] In determining the flow profiles based on data obtained during the shut-in 

period, the well model utilizes thermal characteristics that occur during shut-in. In 
Figures 1 1 and 12, graphs of temperature plotted against distance along the horizontal 
section of the wellbore are provided for the scenarios described above with reference to 
Figures 6 and 7. However, the data graphically illustrated in Figures 1 1 and 12 
represents temperature profiles taken during a shut-in period following the injection 
period illustrated graphically for a first scenario in Figures 6 and 8 and for a second 
scenario in Figures 7 and 9, respectively. The graphs of Figures 1 1 and 12 illustrate 
temperature changes along the horizontal wellbore section at various time points of the 
shut-in. Specifically, the temperature changes are indicated by a graph line 112 
providing a temperature profile at the start of the shut-in period, a graph line 114 
providing a temperature profile at 0.5 days into the shut-in, and a graph line 116 
providing a temperature profile at 1 day into the shut-in period. From the temperature 
profile data, it becomes apparent that the temperature rebounds quickly at intervals of the 
deviated wellbore having a lower rate of injection. On the contrary, intervals with greater 
injectivity, i.e. a greater rate of flow into the reservoir, rebound more slowly. The 
differences in thermal characteristics of the temperature recovery along the length of the 
deviated well enable determination of the injection flow profile. 

[0039] In some applications, the accuracy of the flow profiles can be improved by 

accounting for additional well related parameters. As illustrated in Figure 13, the use of a 
subject model 118 can include additional inputs other then the primary input of 
temperature profiles 120. In this example, the model is utilized or processed on 
automated processor system 62, and a variety of data is fed into the model and processor 
system 62 via, for example, sensors or manual input via input device 68. For example, 
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temperature profile data 120 may be provided by distributed temperature sensor 48. 
Other well related parameters, such as recent history 122, permeability of the reservoir 
124, injection rate 126, injection period 128, and/or thermal conductivity 130, can be 
utilized by model 1 18 on processor system 62 to provide reliable injection flow profiles 
to a well operator. 

[0040] A specific model/algorithm for determining flow profiles based on thermal 

data obtained during injection of fluid may take a variety of physical phenomena into 
account For example, the injection of a cool fluid into a relatively hot reservoir creates 
both a flow of fluid and a flow of heat. Cool or cold fluid moves through the wellbore 
and into the reservoir as heat flows from the reservoir toward the wellbore. A similar 
effect occurs along the wellbore axis in that fluid flows from the heel to the toe, and heat 
flows from the toe to the heel. 

[0041] A practical model to predict the temperature distribution along the 

wellbore when the injection flux is specified, or to estimate the injection flux distribution 
with measured temperature profile can be described as follows. Initial assumptions are 
set forth in the following points: 

(1) The pressure and temperature gradient in the direction parallel to the axis of the 
horizontal well (x) are much smaller than those in the directions (y and z) 
perpendicular to the axis. So, the mass and heat transference along x-direction 
can be neglected. 

(2) The viscosity and density of oil and water are constants. So, the pressure and 
water saturation distribution are independent of temperature. 

(3) Isotropic formation. 

(4) Incompressible injected fluid. 

(5) Open hole completion. 

[0042] The wellbore model in this embodiment further utilizes a wellbore flow 

rate distribution equation and a temperature distribution equation described below. 
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(1) Wellbore flow rate distribution equation 

Let Q^ixyt) denote the injected fluid volume flow rate along the 

wellbore (m 3 /hour), and q^ix^t) denote the fluid volume rate 
injected into the formation for unit length of the wellbore 
(m 2 / hour). The mass conservation of injected fluid in the 
wellbore gives: 

Q,n (*> t) + dx, t) - qwi (x, t)dx = 0 



dx 

Integrating equation (1.1) from the heel to toe yields 



= -^wO>0 (1.1) 



Qmjif) = 2^(0,0 = j gwi (x,t)dx 

0 

QAX'O-QooW-jq^G.tW (1.2) 



0 

where (t) is the total injection rate at the heel of a horizontal well. 





t 

















Wellbore flow rate distribution 



(2) Temperature distribution equation: 
The total heat stored in the tiny element dx is: 

Q well =c,A w T w (x,t)dx 



(1.3) 
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Let ^(x,/) = Itut^t]^—^ denote the heat flow rate from the formation into the unit 

length wellbore and T w (x,t) denote the temperature profile along the wellbore. The 
energy conservation equation for element dx is: 

where q^ and denote the heat rate flowing into and out of the tiny wellbore 

element dx. is composed of two terms: the heat carried by the fluid flowing into the 
element through the wellbore cross-sectional area at x, ^Q^OkjO^VCM)* ^ the heat 
flowing from the formation to the wellbore element through the wellbore surface due to 
heat conduction, q^xj^dx . q^ is also composed of two terms: the heat carried by the 

fluid flowing out of the element through the wellbore cross-section area at x + dx, 
c iQwi i x + dx, t)T w (x + dx, t) , and the heat carried by the fluid flowing out of the element 
through the wellbore surface, c / g w/ (x,0r w (x,/)dx . 



= CiQv* (* + dx, t)T w (x + dx, t) + c tgwl {x, t)T w (x, t)dx . 

Substituting q^ , and equation (1.3) into (1.4) gives: 

d(3*TJ „ (dT\ m , dT w 

- *'. + 2tnpd-z-\ -c l g wl T w =c i A w - 



dx 1 "\ dr w ' " dt 
Substituting equation (1.1) into the above equation yields: 



where c, — heat capacity of injected fluid (J/(mViT)) ; 
— flowing area in the wellbore (m 2 ) ; and 
7] — thermal conductivity of the formation (J/(m. 0 ^ JO) • 
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[0043] Near Wellbore Heat Transportation - Water Injection 

The near wellbore flow regime can be regarded as steady-state radial flow. Consider a 
tiny radial element between r and r + dr . The total heat stored in this element for water 
injection is composed of three terms: the heat stored in the water phase 
Q w = 2m^r<fs w (x,r,t)cjr(x,r,t) , the heat stored in the oil phase 
Q a =2nrdr<j>\L-s w (x,r,t)}: 0 T(x,r,t) , and the heat stored in the rock 
Q r =2wrdr(l-^)c r r(x,r,0 , i.e., 

Q m =Q w +Qo+Q r (L6) 
= Im-ArTix, r, t)[s w (x, r, 0*(c w - c 0 ) + + (1 - <f>)c r ] 

where T(x,r,t) is the temperature distribution in the reservoir, <j> is porosity, c 0 is the 
heat capacity of oil (J/(m 3 :K)) , and s w (x,r,t) is water saturation. 

Let q T (x,r,t) = 1m\ r denote the heat inward radial flow rate, q w (x,r,t) denote the 

dr 

water volume outward radial flow rate, and q 0 (x,r,t) denote the oil volume outward 
radial flow rate. The energy conservation equation for this radial element is: 

where q^ and q^ denote the heat rate flowing into and out of the tiny radial 
element dr . q rin is composed of two terms: the heat carried by the oil and water flowing 
into the element through the inner surface at r, [c w ^ w (x,r,0 + c o^o(^ r ^)l^(^^0» 311(1 
the heat flowing into the element through the outer surface at r + dr due to heat 
conduction, q T (x,r + dr,0 . q^ is also composed of two terms: the heat carried by the 
oil and water flowing out of the radial element through the outer surface at r + dr, 
[c w ^ w (jc,r + dr,0 + c^ o (x,r + dr,0]T(x,r + dr,0 ) ^d the heat flowing out of the tiny 
element through the inner surface at r due to heat conduction, q T (x,r,t) . 

q«n = [c w q w (*>r,t) + c 0 q o (x>r> OFC*^,/) + q T (x,r + dr,0 
S™, = k<7 w (*, r + dr,/) + c 0 q o (x,r 4- dr,t)]T(x,r + dr,*) + q T (r,t) 

dT 

q T {x y rj) = 2m]r—. 
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Substituting the above equations and equation (1.6) into equation (1.7) yields: 
k(W)0(c w _ Co ) + c ^+(i_^) Cr ]^ + ^( Cw - Co )T(x,r,t)^ 



i |~ 2 ar^i 



a* 



a* 



dr 0 dr 

Let (x,t) denote the water front at time t , and choose r^r^ (x,f) . Then: 

s w (x>r,t) = l-s or , q w (x,r,t) = q wi (x,t), q o (x,r,t) = 0 
Near the wellbore, the heat flux also can be considered as a constant, i.e., 

. d( dT\ dq T n 

Equation (1.8) becomes: 

3/ 2nr dr 



(1.8) 



(1.9) 



At the wellbore ( r = r w ), and: 

[d - - O + c> + (i = -^^^[f] . (i-io) 

Combining equation (1.10) with equation (1.5) and choosing c i = c w yields: 



-^=-^0*— ^ . (l.n) 



Let: 



1 



+ f£2- + ( l_0)£r., a 2 =^ 



L^fl + a lw a 2 4? 7 ] 



Equation (1.11) becomes: 



(1.12) 
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[0044] Near Wellbore Heat Transportation - Oil Injection 

For oil injection, the water phase flow rate is g w (x 9 r 9 t) = 0 and water saturation is 

s w (x 9 r 9 f) = . The total heat stored in the radial element is composed of three terms: the 
heat stored in the water phase Q w = lurAnjfs^cJTix.rJ) , the heat stored in the oil phase 
Q 0 =27rrdr0[l-5 wf ]c o r(x,r,O , and the heat stored in the rock 
Q r = 2xrdr(l-(l>)c r T(x,r>t) , i.e., 

Qres=Q*+Qo+Qr 

= 2nr6rT{x 9 r, t)[s^{c„ - c 0 ) + c 0 <f> + (1 - 0)c r ] 

Similarly: 

<lrin = c oQo (*> r > *)T(x, r, 0 + q T (x 9 r + dr, /) 
= Co?* (*, r + dr, t)T(x 9 r + dr,t) + q T (r, f) . 
The energy conservation equation for this radial element is: 
Substituting the above equations into equation (1.7) yields: 

In the steady state flow regime near the wellbore: 

q 0 (x,r,t) = qwl (x,t), — I r— 1 = 0. 
And, equation (1.13) becomes: 

KHc w -c 0 ) + cJ + (l-<l»c r )f = £ ^^f OH) 
At the wellbore ( r = r w ): 

KHc -c.) + cA*(X-*y r f^=-^^{f)^ (1.15) 

Combining equation (1.15) with equation (1.5) and choosing c. = c Q yields: 

l^Bi^-'^'*^-*^^}^- 0 '^ ■ tii6) 
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Lee a l0 =^[l-— |+^+a-^, a 2 = 



1 . 4ot7 



Equation (1.16) becomes: 



[0045] Boundary Condition and Initial Condition 

The boundary condition is specified with the injection temperature at the heel 
(x-0): 

r w (0,0 = 7U(0 (1.18) 
The temperature at the heel 7^(0 can be determined with a wellbore heat 
transmission model, such as the H. J. Ramey model. 
The initial condition is: 

T w (x 9 0) = T R (1.19) 

where T R is the reservoir temperature. 

[0046] Wellbore Temperature Prediction - Forward Problem 

When the reservoir properties (porosity 0 and permeability k), fluid properties 
(density p„andp w , viscosity /u 0 andju w , relative permeability k m andJt^), wellbore 
geometry (wellbore diameter r w , fluid flow area A w , roughness e, and the length of 
perforated section L) are specified, then the injection rate distribution q xvi {x,t) can be 
determined with an analytical model, such as the model established by TUPREP, or a 
numerical model, such as the ECLJPSE100. And thus, with the properly defined 
boundary condition (1.18) and initial equation (1.19), the wellbore temperature profile 
T w ( x >0 can be predicted by solving equation (1.12) for water injection or equation (1.17) 
for oil injection. 
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Let g Q = ^ mj j °^ denote the average injection flux at time t = t 0 , and 



1o 



K c »<lo / 



l + a l0 a 2 



Airr\ 
c Jlo 



Equations (1.12) and (1.17) can be rewritten as: 



dt 



+ a Do (x,t)La, 



dx 



= 0 . 



(1.20) 
(1.21) 



m 1 

The unit of q^ix.t) and g 0 is , so the unit of a w and a n is . 

hour hour 

x 

Let g = — , t D = a J for water injection or t D = a Q t for oil injection. They denote the 
dimensionless variables. The dimensionless forms of equations (1.20) and (1.21 ) are: 

(1.22) 



dt D 



dg 



dg 



Or, both equations can be rewritten as: 



-rf L + «i M fe,f fl )-r^ = 0; (i = o,w) 
ot D dg 

Let g c (t D ) denote the characteristic curve along which the temperature is 

unchanging, i.e., 

d r w =^ + ^^ = o . 

3t D dg e dt D 
Comparing equation (1.24) with (1.25 ) yields: 

dg, 



(1.23) 



(1.24) 



dr. 



(1.25) 



(1.26) 
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Equation (1.26) is the characteristic equation with respect to the partial differential 
equation (1.24). Equation (L26) defines a group of curves, characteristic curves. It can 
be proved that all characteristic curves do not intersect with each other. If one 
characteristic curve crosses the positive q coordinate, then the temperature on this curve 
is specified by the initial condition, i.e., equal to the reservoir temperature T R . Otherwise, 
the curve will cross the positive t D coordinate, and the temperature on this curve is 
specified by the boundary condition T w (t Dp ), where t Dp is the intersection of the 
characteristic curve with the time coordinate. 

[0047] The modeling technique described above enables the determination of 

injection flow profiles based in large part on temperature profiles obtained during 
injection of the fluid. However, the shut-in period also can be modeled such that 
injection flow profiles can be determined based on thermal information obtained during 
the shut-in period. Of course, the data obtained and modeled during the injection period 
and the shut-in period can both be used in determining an injection profile. Furthermore, 
the thermal data obtained when injection is resumed after a shut-in period or the data 
obtained from repeated injection and shut-in periods all can be combined to determine 
and/or verify an injection flow profile. 

[0048] An example of a modeling technique that utilizes thermal data obtained 

during a shut-in period to derive injection flow profiles is described in the following 
paragraphs. First, it should be noted that the temperature profile in the wellbore is 
affected by the fluid convection and the heat conduction between the wellbore and the 
reservoir. Because the thermal behavior of the well depends on the temperature 
distribution around the wellbore, a refined grid block scheme can be used in modeling. 
As illustrated in Figure 14, the model utilizes a grid system 132 that extends in the x, y, 
and z directions. By refining the grid size around the wellbore 30 the temperature profile 
can be stabilized. In other words, the model can utilize a grid system having a grid size 
selected such that further refinement of the individual grid sizes does not affect the 
temperature. 
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[0049] In this type of model, it can be assumed that the temperature distribution 

in the reservoir at the shut-in has the shape of two distinctive regions, one with average 
reservoir temperature and one with the temperature of the wellbore at the shut-in. The 
temperature behavior at the wellbore can be expressed as: 

r »=^fc~<# f - «- 27 > 

The solution of this equation is illustrated in Figure 15 in terms of the dimensionless 
temperature and dimensionless time. Assuming that the temperature front is the same as 
the saturation front, this solution can be used to relate the total injection liquid, water in 
this example, into each interval and the radius of the temperature front is as follows: 

*^=^(l-Sj, (1.28) 

Or, in terms of total injection rate: 

Qtnj^A^Litjx , (1.29) 

This equation is used to estimate the flow profile along the wellbore with t Di determined 
from Figure 15. This method relates injection rate to the temperature in an exponential 
relationship. 

[0050] Thus, a procedure for estimating the injection profile based on thermal 

data obtained during a shut-in period can be summarized as set forth in the flowchart of 
Figure 16. First, the initial temperature of the reservoir, T r , is obtained, as illustrated by 
block 134. The dimensionless temperature for each interval of grid system 132 is 
determined from the temperature log and the reservoir temperature: Td=={T^T s )/(T^TujJ) 9 
as illustrated by block 136. For each unit depth of grid system 132, the dimensionless 
time to can be determined from the graph illustrated in Figure 15, as set forth in block 
138. The dimensionless time for each interval along wellbore 30 is used in the equation 
for estimating the correction coefficient A, as illustrated by block 140. The flow rate into 
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the reservoir for this interval is: A*t Di , as illustrated by block 142, and these flow rates 
can be combined to determine the overall injection flow profile. 

[0051] Accordingly, models, such as those described above, can be used to enable 

the determination of injection flow profiles in deviated wells, such as horizontal wells. 
The use of temperature sensing systems, such as distributed temperature sensors, further 
enable the desired collection of thermal data utilized by the models in deriving accurate 
injection flow profiles. 

[0052] Although only a few embodiments of the present invention have been 

described in detail above, those of ordinary skill in the art will readily appreciate that 
many modifications are possible without materially departing from the teachings of this 
invention. Accordingly, such modifications are intended to be included within the scope 
of this invention as defined in the claims. 
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